Phonon runaway in nanotube quantum dots 
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We explore electronic transport in a nanotube quantum dot strongly coupled with vibrations and 
weakly with leads and the thermal environment. We show that the recent observation of anomalous 
conductance signatures in single-walled carbon nanotube (SWCNT) quantum dots can be 
understood quantitatively in terms of current driven 'hot phonons' that are strongly correlated with 
electrons. Using rate equations in the many-body configuration space for the joint electron-phonon 
distribution, we argue that the variations are indicative of strong electron-phonon coupling requiring 
an analysis beyond the traditional uncorrelated phonon-assisted transport (Tien-Gordon) approach. 
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I. INTRODUCTION 



One of the significant challenges in microelectronics is 
controlling the rapidly increasing thermal budget associ- 
ated with current flow through shrinking devices. Exper- 
imentalQ, S i, 0, H H and theoretical 0, H IH IH 
Il4| investigations are revealing intriguing aspects 
of the mutual effect of electronic and vibronic modes on 
each other. Nanoscalc vibrations tend to couple strongly 
with electronic currents and weakly with their 'macroen- 
vironment', allowing them to be easily driven far from 
equilibrium. Understanding the dynamics of such elec- 
tronically driven phonon runaway processes is crucial to 
the evolution of low-power devices, not to mention the 
novel concepts like molecular motors and phonon lasers 

m 

In this paper, we develop a theoretical treatment of 
current driven nonequilibrium correlated phonon dynam- 
ics in nanoscale systems, and use this approach to analyze 
recent experiments on SWCNT quantum dots (QD) 
Using a rate equation for correlated transport in the full 
many-body eigenspace of the coupled electron-phonon- 
lead-bath system, we explain novel spectroscopic features 
such as the anomalously large absorption sidebands aris- 
ing from phonon runaway in suspended. Coulomb Block- 
aded nanotubcs (Fig. 1). 

Our model also explains semi-quantitative features of 
the experiment such as the amplitude variation of the 
Coulomb Blockaded conductance peaks and their phonon 
sidebands, as a function of injected current. However, 
our model predicts a linear variation in phonon popula- 
tion with current, in contrast with experiments that show 
a quadratic variation P, |2| ■ Wc argue that a possible ori- 
gin of this discrepancy is because our model explicitly 
incorporates the effect of strong electron-phonon correla- 
tion that is characteristic of these experiments, while the 
experimentally extracted variation was accomplished by 
employing a traditional Tien-Gordon analysis that im- 
plicitly treats the phonon contribution only through its 
mean-field oscillating potential acting on the electronic 
subsystem [Tij . 



II. MODEL 

We use a model Hamiltonian Hd for a quantum dot with 
onsite energies e^, Coulomb interaction energy Un', vi- 
bronic modes at energy tiWj and electron-phonon cou- 
pling \ij (Fig. [2). The total Hamiltonian including the 
contacts (He), the phonon bath (Hb) and their couplings 
with the dot [Hoci Hdb) is: 
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Hd + He + Hb + Hdc + Hdb (1) 

Hel + Hph + Hel-ph (2) 

j 

^ Ay?kjj7ij(a] + Qj) 

£kadl^dka (3) 

k.aeL.R 

Y,fiujib\bi (4) 

Tika{d\^Ci + dkacl) (5) 

i,k,a 

J2^Abl + bi)ia]+a,) (6) 



where, (c) and d^ {d) are the electronic creation (de- 
struction) operators for the dot and the leads, = cjci 
and (a) and (b) are the phonon creation (destruc- 
tion) operators for the dot and the phonon bath respec- 
tively, r, ek and k represent respectively the dot-contact 
coupling, the contact bandstructure and the coupling be- 
tween the dot vibrations and the thermal bath). 

The electron-phonon coupling is eliminated using a 
standard unitary polaronic transformation H = e^He~^ 
[iTj . where S = j ^ij (o-] ~ o-j)ni- This transformation 
renormalizes the onsite (e^) and Coulomb {Uu') energies 
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for the dot: 

Hd = ^ Ejcjcj + ^ ^ ?7ij'nj7V + ^ ^ja]aj (7) 
h = e.-^^X^hujj (8) 

3 

Uii' = Uii' — 2 ^ ^ XjjXj'j (9) 
j 

where the 'dressed' electron or polaronic annihilation op- 
erator 



e-^ce"-^ = cX 



(10) 



with X = cxp^^- Xij{aj — aj) denoting the generator of 
the polaronic shift and 



e a,e 



(11) 



while preserving the electronic number operator fii = rii. 
At this point, we will simplify the model by considering 
only a single phonon mode, denoted by the index 1. 

Current flow in this system involves single elec- 
tron transitions between many-body states |e]yg,fc) (fc 
phonons and ith electronic level in the iVe electronic sub- 
space) of the quantum dot, with matrix elements for the 
electronic destruction operator: 



e->^"^'\~Kif-^ ^Cl-P{^l) for k>p 
e-^-/2(-A,i)f-\/t£P-^-(A2j for p > k 



(12) 



where, Ci\ej^^,k >~ \e 



Nf,_i,k > and are the associ- 
ated Laguerre polynomials. For tunneling rates Ta/h (a 
= L(left)/R(right)) (Fig.O, the transition rates are: 



fj£{\e%„p))-£{\e'-j,,_,,k))-r,\q\Vappi (13) 
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FIG. 1: (a): Schematic of STM measurement on CNT QD; (b) 
and (c): Experimental observation and theoretical calculation 
for CNT QD at lower and higher current level {solid line: 
conductance, dashed line: (Nph) in (c)). 1,2,..., 6 denotes the 
main Coulomb peaks. The parameters for the calculation are- 
ei = £2 = 24 meV, h = h = h = h = 54 meV, U = 27 meV, 
Ef = 0, huji = 11.5 meV, An = A21 = ■■ ■ = Aei = 1.6, 
Tl = 7.5 X 10"'' eV, /3 = 5 X 10"® eV, = 0.6 and T = 5K. 



/i.fl, are the contact Fermi functions with electrochemi- 
cal potentials = Ep, = Ep — \q\Vappi at tempera- 
tures Ti,i{, while the electrostatic voltage division factor 
77 represents the fraction of the applied bias acting on the 
levels. Here £{\)) stands for the energy of the eigenstate 
I) of the dot. 

In our analysis, electrons enter or leave the dot by emit- 
ting or absorbing phonons, which in turn are coupled 
with their thermal environment, held at a bath tempera- 
ture Tb, with an escape rate P/h (Fig.|21), that maintains 
a Boltzmann ratio between emission and absorption pro- 
cesses in the contact 
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/3, 



(14) 



The state transitions in the joint electron-phonon 
many-body space are described by the master equation 



,k) 



dt 



r,N' 



(15) 



together with the normalization condition for P|es j,^, 
where 7^ = + + 7^P'*. 

Solving the rate equations at steady state gives us the 
current / and the steady state population of dot electrons 
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FIG. 2; The dot is electrically connected to the left (right) 
contact (with electron tunneling rates Tl,r/K) and mechani- 
cally to the phonon bath (with a phonon escape rate of P/h). 
The dot has electronic degrees of freedom ei and phonon de- 
grees of freedom hujj = 1, 2, 3, . . .) with coupling Xij. 



(Nei) and phonons (Nph): 

N^,N'^,r,s,k,n 

(Npk) - ^^'^■ive.fc) (16) 

where, sgn{x) is the signum function. This also gives us 
the phonon generation rate Q°' by the current and the 
phonon extraction rate X^^ (see appendix) by the bath: 



s,N^,k r,N^,n 
s.N^.k r.N'.n 



ph 



(17) 



where, iVp' is the equilibrium phonon occupancy at the 
bath temperature. 

The rate equations are schematically explained in 
Fig. 12 The inputs to these equations are the electron 
tunneling rates TL,R^ the phonon escape rate /?, the elec- 
tronic energy configuration e^, the phonon energy fiWj, 
the charging energy U and the electron-phonon coupling 
Ay . In our analysis we will set the lead and bath tem- 
peratures T]^ ji and Tg to be equal to the ambient tem- 
perature. 



III. NON-EQUILIBRIUM PHONON 
OCCUPATION IN CNT QD 

A. Peaks in the conductance spectrum due to 
phonon-assisted tunneling 

We will now apply our many-body rate equations to an- 
alyze recent experiments Q on phonon-assisted tun- 
neling in suspended SWCNT quantum dots. The experi- 
ment shows several striking features such as anomalously 
large phonon absorption peaks at low temperature (5 K) 
and a monotonic increase in phonon sideband amplitude 
with current (Fig.^b and c). Conductance peaks are ob- 
served in groups of four, suggesting consecutive doubly- 
spin-degenerate electronic levels in the Coulomb Block- 
ade regime. Sidebands are attributed to phonon-assisted 
tunneling through the radial breathing mode (RBM) of 
the SWCNT, known to predominate at low bias through 
its effect at the bottom of the conduction band ■ 

We will model the dot with a single phonon mode {hu\) 
corresponding to the RBM and six electronic energy lev- 
els, - two lower energy degenerate levels (ei and I2) and 
four higher energy degenerate ones (53, £4, £5 and eg). 
The occurrence of the singly-degenerate and the doubly- 
degenerate discrete electronic levels in CNT QD has been 
observed earlier |23| . We will assume that the tun- 
neling rate Ta/h and the phonon escape rate f3/h are 
dispersionless, and use the same value for the couplings 
(All = A21 = • • • = Agi) of the RBM to all the electronic 
levels [13. 

In order to match experimental results, we find it nec- 
essary to consider electron addition levels only, and not 
electron removal levels lying below the equilibrium Fermi 
energy 0, 0. An analytical estimate of the height of 
the first two conductance peaks, corresponding to two 
degenerate levels at 61^2, explains the justification be- 
hind this assumption. It can be shown that the height 
of the first two electron removal peaks should be propor- 
tional to 2rirfi/(2rR + r^) and 2rLrR/{rB. + Tl) - 

2rLrfl/(2rfl + Fl) [23, assuming that the left contact 
injects electron, the temperature is very low, and ignor- 
ing phonon sidebands, - assumptions which are consistent 
with experimental conditions. However, this predicts, for 
F/i ^ Fi, that the height of the second peak should be- 
come zero in contrast with the experiment (Fig.^). The 
height of the first two electron addition peaks, under the 
same assumption, can be shown to be proportional to 
2FLFfl/(Ffl+2FL) and 2FiFfl/(Ffl,+FL)-2FLFfl/(Ffl-t- 
2Fi), which predicts, for F/? ^ F^, that the peak heights 
should be equal, in agreement with experiment (Fig.^). 
This suggests that the observed peaks arise from electron 
addition rather than removal. 

Next, let us estimate the input parameters. The 
phonon energy huji was measured to be 11.5 meV (see 
which together with the separation between the 
main Coulomb peak and its first phonon emission side- 
band, yields a voltage-division factor 77 « 0.6. The po- 
laron renormalized charging parameter U is estimated 
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FIG. 3: Variation of conductance with the current at the 
second coulomb peak (at Vappi ~ —0.1 V) and it's associated 
first phonon sidepeak (at Vappi = —0.125 V). inset: Variation 
of the no. of phonons with current at the above mentioned 
conductance peaks. The parameters used to generate these 
results are; ii = €2 = e — 8 meV, U = 12 meV, Ef = 0, 
huji = 12.5 meV, An = A21 = 2.1, Tl = 5 x 10"^ eV, Fr = 
0.00005rL to 0.5rL, /? = 5 x lO^"' eV, r; = 0.2 and T = 5K. 
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FIG. 4: Simple one single electron level model do not match 
the experimental result at all. The values of the parameters 
used are: ei = 0.05 eV, ftoji = 12.5 meV, An = 1.7, Fl = 
1 X 10-^ eV, Vr = 0.00005Fl to 0.5Fl, /3 = 5 x 10"^° eV, 
77 = 0.5 and T = 5 K 



to be 30 meV, as extracted from the separation between 
consecutive Coulomb peaks originating from the same de- 
generate set of levels. From the estimate of the phonon 
decay rate and the Q factor, f3 was determined to be 
~ 10^® eV. In our calculations we varied the tunnel- 
ing rates V ^ and clcctron-phonon coupling constant A^i 
(i = 1, 2, . . . , 6) to match the experimental conductance 
levels. 

Figs.^b and c show a comparison between the exper- 
iment and our calculations with the above parameters. 
Six levels seem sufficient to capture the essential physics, 
including the increase in number of phonons between 
emission and absorption sidepeaks arising from the cor- 
responding increase or decrease of phonon occupation at 
those bias points (Fig.^c). The calculated phonon num- 
ber significantly exceeds the equilibrium value after each 
emission peak and drops considerably after absorption, 
indicating strongly correlated nonequilibrium phonon dy- 
namics in this system. The weak phonon-substrate cou- 
pling P for suspended tubes leads to a phonon bottle- 
neck whereby the current emits phonons faster than they 
are conducted away, leading to anomalous low tempera- 
ture absorption peaks that even exceed the corresponding 
emission peak heights on occasion. 



B. Variation of conductance with current 

Experiment Q shows a prominent quadratic dependence 
of the phonon occupancy on current. This variation is ex- 
tracted indirectly from the observed variation in height of 
a main conduction peak and it's associated first phonon 



sidepeak with current (Fig. Ql), by applying a traditional 
Tien-Gordon model to fit it Note that these re- 

ported experimental results are from a different CNT QD 
sample from the one corresponding to Figs.^b and c. So 
we will use a model with a different configuration of pa- 
rameters to capture the essential characteristics of the 
experimental observations. In order to ascertain the con- 
ductance variation with current, we have replotted the 
reported normalized conductance, [dl / dVappi) / {I /Vappi) 
vs. current, / to conductance G = dl / dVappi vs. current 
/, using the known bias values Vappi (-0.1 V at main peak 
and -0.125 V at first sidepeak) and current values / at 
those peaks obtained from 0. The experimental obser- 
vation shows that the main conductance peak initially 
increases with current and thereafter tends to saturate. 
On the other hand, the conductance at the first side peak 
gradually increases with current exponentially within the 
entire current range of the experiment. 

To explore this variation we further simplify our model 
to just two degenerate single electronic levels at energy 
e. The simulation produces two Coulomb peaks and 
their associated phonon sidepeaks. To mimic the exper- 
imental procedure, we increased the current by varying 
the tunneling rate of the STM tip. The calculated 
current dependences of the conductance at the second 
main peak (due to the direct tunneling at energy e + ?7) 
and its associated first phonon sidepeak at higher en- 
ergy [i -\- U -\- HuJi) agree with the experimental trend 
(Fig. 13). However, we do not find (Nph) varying as P 
at the conduction peaks {Fig. ^ inset), but instead al- 
most linearly within the experimental current range. The 
disagreement between our result and the claim made in 
111 stems from the fundamental difference between our 
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FIG. 5: Variation of conductance with current for current 
values higher than the experimental range, inset: The vari- 
ation of conductance at the coulomb peak in the absence of 
electron phonon coupling. The values of the parameters used 
are the same as those of Fig. |21 (except that their is no cou- 
pling between electrons and phonons for the results in inset). 
The region inside the dotted box is the range of experimental 
observation in Fig. |3] 



model and the Tien-Gordon model. Unlike our model, 
the Tien-Gordon model neglects the correlation between 
the electronic degrees of freedom and the phonon de- 
grees of freedom by replacing the electron-phonon cou- 
pling Hamiltonian Ifoi-ph with a mean-field oscillating 
potential V = {5Hc\-p\xl 5n) which is proportional to the 
average position operator. For a particular phonon mode 
of frequency lo and in the absence of leads or baths, the 
position operator and thus the averaged phonon poten- 
tial is simply proportional to coswt. Since the Bohr fre- 
quency of the electron in turn depends on this oscillat- 
ing Hamiltonian, the resulting time-evolution operator 
and thus the electronic spectral functions end up with 
independent Fourier components whose spectral weights 
depend on Bessel functions, in other words, the Tien- 
Gordon expression. Crucial to the derivation of the Tien- 
Gordon model is the replacement of the electron-phonon 
coupling Hamiltonian by its mean-field time-dependent 
component ignoring leads and baths, which eliminates 
all correlation effects. This assumption is clearly in- 
consistent with our model that shows strongly corre- 
lated phonon dynamics, captured by the dimensionless 
electron-phonon coupling constant A that is typically 
greater than unity, as mentioned in our figure captions. 

It is worth explaining why a single level does not suf- 
fice to capture this current variation. We find that the 
conductance of the main Coulomb peak, irrespective of 
the parameters chosen, varies linearly with current and 
does not saturate. Since the current prior to the main 
Coulomb peak is zero, the height of the main conductance 
peak ~ Icoui/'^kBT , where Icoui is the current plateau im- 
mediately past this voltage, the main broadening coming 
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FIG. 6: Variation of conductance with current for four de- 
generate single electron levels at the four Coulomb peaks and 
their associated first phonon sidepeaks. The values of the pa- 
rameters used are: ei = £2 = es = £4 = e = 20 meV, ?7 = 30 
meV, Ef — 0, huji = 12.5 meV, An — A21 — X31 — Xu = 2.1, 
Fi = 5 X 10"'' eV, Tr = O.OOOOSFl to 0.5Fi, /? = 5 x 10~^° 
eV, 77 = 0.5 and T = 5K. 



from temperature. This result, consistent with a test sim- 
ulation with one single electron level (Fig0J, also persists 
for the first in a series of Coulomb peaks in a multilevel 
system for the same reason. In all these examples, the 
main conductance peak can only vary linearly with the 
current, in contradiction with experiment, prompting us 
to look at higher peaks in a more complex, multileveled 
system. This requires us to adopt a more complex model 
with two single electron levels and to look at the varia- 
tion at second Coulomb peak and it's associated phonon 
sidepeak. The results reported in the experiment (Figs. 
4a and 4b of ) were later confirmed to involve the vari- 
ation of second Coulomb peak and its associated phonon 
sidepeak pl |. 

We used our doubly degenerate single electron level 
model to explore the variation of the conductance for 
current values beyond the experimental current range, 
to project the characteristics under higher STM set cur- 
rent. The results are shown in Fig.|Sl The conductance 
values at both the Coulomb peak and the first sidepeak 
show a non-trivial variation with current. The variation 
of the conductance at the Coulomb peak is compared 
to that in (Fig. |31 inset) absence of electron-phonon cou- 
pling. We find that the coupling of electronic and phonon 
degrees of freedom significantly affects the Coulomb peak 
characteristics. In particular, while the trends match the 
experimental data up to the experimental current lev- 
els, we notice that at higher currents the main peak goes 
through a kink followed by a rise, while the sidepeak am- 
plitude tends to saturate. 

Finally we look at the conductance variation with 
current for a quadruply degenerate single electron level 
model for a CNT QD that occasionally shows doubly 



spin degenerate shells. The conductance variation at 
each Coulomb peak and its associated first higher en- 
ergy phonon sidepeak are shown in Fig. The first 
Coulomb peak varies almost linearly with current (sim- 
ilar to Fig. 01), as we argued earlier. The conductance 
variations at the other three peaks show the same qual- 
itative trends as Fig. |31 From the qualitative similarity 
between Figures 0] |S1 and El we sec that the qualitative 
features of these variations are robust with respect to the 
values of the parameters that went into the model. 



IV. CONCLUSION 

We have studied strong electron-phonon coupled dy- 
namics under nonequilibrium conditions using a many- 
electron rate equation approach, focusing on recent ex- 
periments on suspended nanotubes. From our calcu- 
lations we argue that phonons in weakly coupled sys- 
tems can be readily driven far from equilibrium, lead- 
ing to anomalous signatures in the corresponding con- 
ductance spectrum and temperature dependences. In- 
teresting extensions of this work would involve studying 
the imputations of nonequilibrium phonon dynamics for 
energy dissipation in nanoscale systems, and of main- 
taining phonon coherence by tuning their decay rates 
through their mechanical coupling with the substrate, 
and through phonon-phonon interactions controlled by 
the inherent nonlinearity of the lattice. 
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VI. APPENDIX 



When A = Ci we have: 

[S,Ci] = ^Xhj{a^j - aj)[clch,Ci\ 



= ^ (a] - aj)[clchCi ~ Ciclch] 



h,3 



hChl 



h,3 



where 5 is the kroncckcr delta function. Then, 
[S, [S, c^\] = [S, X Xrj [a] - aj)ci] 

= Yl ("1 ~ "j )] H [^hj{a] 

,3 

X Ay (aj - a]) 



(19) 



[S, [S,c,]] 



2 



[(ich,C^] 



(20) 



and so on. Putting all these results in equation 1181 we 
finally get eq. ^1 When A = Uj we have: 

[S,aj] = y^^niXihiaj ~ ah,aj] 

i.h 

= '^niXihiahttj - a\aj - ajUh + a^al) 

i.h 
i.h 

[S^aj] = J2''A^J (21) 



So, afterwards [S,[S,aj]] = [S,[S,[S,aj]]] = • • • = 0. 
Putting all these commutators in equation ^1 we get 
eq.im 



In this section we show the derivation of some impor- 
tant analytical results. The derivations i nlVI Al and lVI Bl 
which are already dicussed in ref. [l^ . Il7j. have been 
worked out here to maintain continuity. 



A. Derivation of Equation 1101 and 1111 

The terms like e^Ae~^ are evaluated using: 

e^Ae-^ ^A+[S,A] + (1/2!) [5, [S, A]] + ... (18) 



B. Derivation of Equation 7, 8 and 9 

Applying the polaronic transformation: 

Hd = Hel + Hph + Hel-ph (22) 

Now the electronic portion: 

Hel = '^<^ini + - X Uwriifii' (23) 

i i.i' 
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The phonon portion: 
Hph = e HphC 



By simple algebra: 



j i i' 

hujja^^aj + hujjXijXiijn.iriii 
hojjXijTiiaj — huijXi'jniiOj 



E 

and similarly 



(A.i)' 



/fc! 



1/2 



|e^,,fc-Z) (27) 



\lj ILi U>j 
1,] t'.J 

Hph — hujjo},^aj + hhjjXijXiijh 



P 



— m) ! 



1/2 



(28) 



E A^^j^j - E ^jA,jnj(a] + a^) (24) p^.^^ equations EHlEZIand [23 



where we have used riirii = Ui = fii. Now, finally, the 
electron phonon coupling term: 

Hel-ph = Hel-phG~^ 

Hei-ph = Xjjhw jni{a\ + g.,-) 

= E f^]Xijni{a}. + - E Ai'j'^I' ^ E Ai'j'^i') 
= E f^jXtjni{a}- + Oj) - 2 E f^jXijXi, jfiifiii 

Hei^ph = E Ajjni(a] - 2 E Xjjhwjh^ 

(25) 



(e^^,_i,p| £, |e^e,fc) = e-'-/^EE 

rn=0 ;=0 

(p!fc!)i/2 



(A,i)'(-A,i)^ 



m\l\ 



^p—rii.k—l 



(29) 



For the case of fc > p, to remove the kroncckcr delta 
function, we substitute I = k — p + m in Eq. 1291 

m\[k-p + m)\ 

(p!fc!)i/2 



1,3 *J 

2 E hwjXijXiijfiifii' 



{p — my. 

(-l)"'(A.i)'-^(Afi)'»(fc!p!)V^ 



Substituting equations (23 ISl f^i^d I^S] in equation [23 we 
get the results of equations 7, 8 and 9. 

C. Derivation of Eauation ll2l 

We will start from: 

Now since Ci\e%^,p >= |e^g_]^,p > it follows that: 

(e^e_i,p|cj;|e^g, fc) = (e^g,p|X|e^g, fc) 
For one phonon mode X = exp[Aii(ai — a\)]. Then: 

Using e-'*'i(°i~°i) = g-AiiaJgAiiaig-A^i/2. 
( ^Afe-liPlCil^Afei ^} — 

e-^-/2(e^„p|e-^-'^Je^-'^Me^„fc) 

(26) 



E 

m— 

e-^-/^(A.a 



?7!,!(fc — p + m)!(p — m)! 



^ ^ (-l)"'{(fc-p)+p}!(Afi)- 
m!{(fc-p) + m}!(p-m)! 

The summed series is nothing but associated Laguerre 
polynomial Cp^P{Xfi). So, finally we get: 

( eNe-i^P\ct\e%e,k) = 

e-^?i/2(Aa)'^-^'V^4-^(A|i) for fc > p 

Following the exact same procedure we can show that, 

( 6jVe-l;P|Ci|eAre, fc) = 



e-^-/^(Aa) 



"-'■y^^r'(A?i) iorp>k 



D. Derivation of Equation 1171 

From the equations 1141 we can see that the phonon 
bath induces a transition between only those two states 
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which have the exact same electronic configuration and 
phonon number differing by 1. So we can write: 

s,Ne,k r.Ne' ,n 



s,N^ k=l 



p 



s,N^ k=0 



ph 



fc - (fc + l)cxp[- 



lEE^i 



So, from Eq. 1161 and the normalization condition: 



^ph 



(Nph) ( 1 



iv7 



cxp[ 



ph 



f3\ {n,,)^n;^ 



ph 



ph 



where, = [exp[ ^^^/(fcBT)] - 
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